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ABSTRACT 

The gravitationally lensed quasar HE 1 104-1805 has been observed at a variety of wavelengths ranging from 
the mid-infrared to X-ray for nearly 20 years. We combine flux ratios from the literature, including recent 
Chandra data, with new observations from the SMARTS telescope and HST, and use them to investigate the 
spatial structure of the central regions using a Bayesian Monte Carlo analysis of the microlensing variability. 
The wide wavelength coverage allows us to constrain not only the accretion disk half-light radius r\n, but the 
power-law slope £ of the size-wavelength relation ryn ex As With a logarithmic prior on the source size, the 
(observed-frame) /?-band half-light radius log(r 1; /2/cm) is 16. 0!^, and the slope £ is 1.0S- We put upper 
limits on the source size in soft (0.4-1. 2keV) and hard (1.2-8 keV) X-ray bands, finding 95% upper limits 
on log(r 1/ / 2 /cm) of 15.33 in both bands. A linear prior yields somewhat larger sizes, particularly in the X- 
ray bands. For comparison, the gravitational radius, using a black hole mass estimated using the H/3 line, is 
log(r ? /cm) = 13.94. We find that the accretion disk is probably close to face-on, with cos/ =1.0 being four 
times more likely than cos/ = 0.5. We also find probability distributions for the mean mass of the stars in the 
foreground lensing galaxy, the direction of the transverse peculiar velocity of the lens, and the position angle 
of the projected accretion disk's major axis (if not face-on). 

Subject headings: accretion, accretion disks — gravitational lensing: micro — quasars: individual 
(HE 1104-1805) 



1. INTRODUCTION 

The detailed structure of the innermost regions of active 
galactic nuclei (AGNs), between 1 and ^1000 gravitational 
radii from the central black hole, has remained observation- 
ally elusive. The gravitational radius of a black hole of mass 
10 9 M Q is about 10 AU, so these scales cannot even remotely 
be resolved at cosmological distances. The ultraviolet (UV), 
optical, and near-infrared continuum is thought to com e from 
a geometrically thin, optically thick accretio n disk ( |SIiakura| 
|& Sunyaev|[T973| [Novikov & Thorne|[T973) . Although this 
model has been quite successful in explaining the X-ray spec- 



tra of stellar-mass black hole binaries (e.g., McClintock et al. 

I20TT), 

WKb 



the same has not been true of the UV/optical spectra 
Ns, in p art because o f complications arising from line 
emission (e.g., Blaes et al. 2001). The X-ray continuum is 
non-thermal and is thought to arise from the inverse Compton 
scattering of disk photons by a corona of hot electrons (e.g., 
Reynolds & Nowak 2003). The spatial structure of the X-ray 



corona is not known. In addition, in many AGNs the presence 
of the iron Ka emission line indicates t hat X-rays are being 
reflected from the accretion disk (e.g., [Fabian et al.|1989||Laor| 
TWT) . 

The gravitational microlensing of lensed quasars has proven 
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to be an effective tool for measuring the properties of quasar 
accretion disks, and is starting to become useful for studying 
the X-ray corona as well. The time-dependent microlensing 
magnification (or demagnification) of one or more images of 
a lensed quasar is moderated by the finite size of the source, 
which smooths the complicated caustic pattern of microlens- 
ing magnifications as the quasar passes over it. This allows us 
to use the microlensing magnifications to estimate the source 
size, and such work has shown that in general the accretion 
disks are larger than would be expected fro m either thin disk 
modeling or tota l flux argument s (Pooley et al. 200 7) |An 
guita et al|[2008l IMorgan et al. l [2~010||Hainlrne et al.||20l' 
Jimenez-Vicente et al. 2012l Since the effective tempera- 
ture of the disk depends on radius, its apparent size depends 
on wavelength, leading to a chromatic dependence of the mi- 
crolensing magnification, with the same quasar image experi- 
encing larger variability at blue wavelengths than at red wave- 
lengths. Several studies have used this to constrain the power- 
law slope of the size-wavelength relation, and the results have 
been consistent with each other and with the thin disk predic- 
tion that the size goes like the 4/3 power of the wavelength, 



mostly because of their large uncertainties (Poindexter et al. 
2008, Bate et al. 2008; Eigenbrod et al. 2008, Floy d et al. 
2009 1 Blackburne et al. 2011b; Mosquera et a l.|2011| l. FinailyT 



efforts to put upper limits on the size of the X-ray regions have 
xessful jPooley et al |2006] [20071 jChartas et al. 
al. 2010), and recently there have been attempts 



also been successful 

2009; Dare? 

to constrain the direct and reflected components' sizes sepa 

rately using color cuts or spectral decomposition (|Chen et al. 



201 1 , Blackburne et al. 201 la; Morgan et al. 2012 Chen et al. 



2012b; Chartas et al. 2012; Mosquera et al. 2013). 

In this paper we use infrared (IR), optical, UV, and X- 
ray photometry of the lensed zs = 2.32 quasar HE 1 104-1805 
(Wisotzki et al.|[T993] |1995| l to quantitatively constrain the 
properties of its centralemission regions. HE 1104 is lensed 
by a foreground early-type galaxy at redshift z L = 0.73 into a 
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pairofimages separated by 3 '.'2, The mass ofthe central black 
hole in HE 1104 has been investigated using the widths of 
the emission lines Civ, H/3, and Ha, yielding mass estimates 
log(M B H/ M Q ) = 9.37±0.33 8.77±0.30,and9.05±0.23,re- 
spectively ( |Peng et al 120061 |Greene et a"T1|2010[ |Assef et al.| 

h7 



TABLE 1 
SMARTS Light Curves 



201 1) . We adopt the R(j mass ol'|Assef et al.|< |2011| fortius 
paper. We compare the light curves to microlensing simula- 



tions using the Bayesian Monte Carlo method of Kochanek 
(2004) and |Poindexter & Kochanek| ( ]2010"a|b"] l, which allows 
us to derive posterior probability distributions for our param- 
eters of interest, which include the half-light radius of the ac- 
cretion disk and the X-ray emission regions, the slope with 
which the radius changes with wavelength, and the inclina- 
tion of the disk, as well as the mean mass of the stars in the 
foreground galaxy. 

2. MULTI-WAVELENGTH DATA 

We combine data from the literature with new photometry 
from the Small and Moderate Aperture Research Telescope 
System (SMARTS) and Hubble Space Telescope (HST) to 
create a lightcurve spanning about 19 years and from the near- 
IR to X-rays in wavelength. 

The previously-published data in our light curve come from 



Courbin et all ( [T998] J, K),\Remy et al.](fT998] V, I, g),|Lehar 
et al. (2000 , H), |Gil-Merino et aT[ (|2002, B), Schecht er et al 



V), Ofek&Maoz (2003, 



RAC 3.6 pm), and Munoz 



(2003), Wyrzykow ski et al.| ( |200 3 
R), Poindexter et al. (2007p/, K 

|etal.| ( T2011[ F330W, F435W, V, F625W, /). Where applica- 
ble, we have used the shorthand V, I, or H for the HS T filters 
F555W F814W , or F160W, respectively. Although Poindex- 
|ter et al.| ( |2007| l report flux ratios from the Spitzer Space Tele- 
scope at several mid-IR wavelengths, we only use the 3.6 /.im 
ratio, as it is the most likely to originate in the accretion disk 
rather than a dusty torus. The mid-IR flux ratios are all nearly 
identical, indicating that the source is large enough at all 
these wavelengths for microlensing not to be important. We 
also note that we have actuall y taken the HS T F160W mag- 
nitudes attributed to Poindext er et al.| ((2007) from the CAS- 
TLES database, since due to a typographical error this paper 
reports the values of Lehar et al. (2000) rather than its own 
measurements. We have confirmed that the CASTLES magni- 
tudes correspond to the same observations (E. Falco, private 
communication). 

2.1. Optical data 

In addition, we use data from 8 seasons of monitoring by 
the ANDICAM camera (DePo y et al |2003[ l on the SMARTS 
telescope, primarily in the R and J bands, with some data in B, 
V, and /. These images are bias-corrected and flat-fielded by 
an automated pipeline, and we stack the three to six images 
obtained on each night of observation to improve the signal to 
noise ratio and reject cosmic rays. We reject epochs with bad 
seeing (FWHM > 2."0) or high sky levels indicative of clouds 
or excessive moonlight. We measure the fluxes of the quasar 
images using the point spread function (PSF) fitting method 
described by Kochanek et al. (2006). The first 3 seasons of th e 
R and J light curves are reported by Poinde xter et al.|f2007| l, 
but for convenience we report them in their entirety, together 
with the B, V, I, and J data, in Table [1] Since we have re- 
analyzed theimages^someof the magnitudes differ slightly 
from the Poinde xter et al.| ( [2007[ ) values. 

HE 1104 has a fairly long lensing time delay between its 
two images, which means that we must be wary of intrinsic 



HJD- 2450000 


A 


B 


Filter 


2976.806 


+0.11 ±0.01 


+1.67 ±0.03 


J 


2976.807 


-0.01 ±0.01 


+1.38±0.02 


R 


2985.815 


+0.13 ±0.01 


+1.56 ±0.03 


.1 


2985.818 


-0.02 ±0.01 


+1.39±0.01 


R 


2993.822 


+0.11 ±0.01 


+1.53 ±0.02 


.1 


2993.826 


-0.03 ±0.02 


+1.40 ±0.02 


R 


3000.816 


+0.10±0.01 


+1.63 ±0.02 


.1 


3000.820 


-0.05 ± 0.02 


+1.40 ±0.02 


R 


3009.812 


+0.06 ±0.01 


+1.61 ±0.02 


.1 


3009.815 


-0.09 ±0.01 


+1.41 ±0.02 


R 



NOTE. — Light curves are in uncalibrated magnitudes. 
This table is published in its entirety online. A portion is 
shown here for guidance regarding its form and content. 



quasar variability conspiring with this delay to mimic mi- 
crolensing variability. F or this work, we adopt a delay of 
162.2 days ([Morgan et al.|2008|l, 



~aLl2 
mely 



with image A leading. For 
datasets (specified uniquely by publication, observatory, and 
filter) with many observations, we use linear interpolation on 
the light curve of image B to remove the time delay. We do not 
extrapolate outside the bounds of any dataset, and we limit the 
interpolation to dates within 30 days of an actual data point, 
in order to avoid interpolating across seasonal gaps. We also 
avoid comparing light curves from different datasets, in or- 
der to avoid systematic calibration errors. Because of this, we 
do not have to worry about the flux calibration of our data, 
working exclusively with magnitude differences, or flux ra- 
tios between the two quasar images. After interpolating, we 
bin sequential pairs of observations if they are separated by 
30 days or less. This is just to keep our light curve short, 
since our simulation software's memory requirements grow 
linearly with the number of epochs. Several datasets, particu- 
larly the early ones, have single-epoch observations, so we are 
forced to use the original (not time delay-corrected) flux ra- 
tios, and to add to their uncertainties some estimate of the er- 
ror caused by intrinsic variability. For this estimate we use the 
root-mean-square (RMS) ofthe difference between the delay- 
corrected and uncorrected /?-band light curves of image B. 
This value comes out to 0.1 1 mags, so for the affected obser- 
vations we add 0.078 mags in quadrature to the uncertainties 
for images A and B, dividing the total between the two. 

Several datasets, particularly the earlier ones, are too 
sparsely sampled to allow linear interpolation according to 
our standard procedure, but nevertheless have better cover- 
age than the single-epoch datasets. In these cases we choose 
pairs of epochs that are fairly well-matched after correcting 
for the time delay, and we associate image A from one epoch 
with image B from another, augmenting their error bars as in 
the single-epoch cases, but by only half the amount. In par- 
ticular, we apply this method to a previously unpublished pair 
of observations of HE 1 104 in the B band from the Instituto 
de Astrofisica de Canarias (IAC) 80-cm telescope at the Ob- 
servatorio del Teide. The respective instrumental magnitudes 
of images A and B were -1.179 ±0.016 and 0.509 ±0.040 
on 19 January 1999, and -1.261 ±0.036 and 0.445 ±0.100 
on 18 March 1999. We construct a single magnitude differ- 
ence from these two observations, adopting the earlier value 
for image B and the later value for image A. 



We give the V-band dat a of Schechter et al. (2003) and 
Wyrzykowski et al. ( |2003[ ) special treatment in two ways. 
First, due to its dense sampling we bin it four observations 
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TABLE 2 
HST F275W LIGHT CURVES 



HJD- 2450000 


A 


B 


Star a 


Starb 


Stare 


5169.817 


17.529 ±0.002 


18.085 ±0.002 


16.094 ±0.001 


16.014± 0.001 


18.170±0.002 


5227.008 


17.184±0.001 


18.123 ±0.002 


16.074 ±0.001 


16.004 ±0.001 


18.637 ±0.003 


5270.614 


17.113±0.001 


18.219 ±0.002 






18.831 ±0.003 


5340.061 


17. 126± 0.001 


18.589 ±0.003 


16.113±0.001 


16.043 ±0.001 


18.344 ±0.002 


5382.058 


16.904 ± 0.001 


18.327 ±0.002 


16.086 ±0.001 


16.008 ±0.001 


18.223 ±0.002 


5505.015 


17.348 ±0.001 


18.328 ±0.002 


16.094 ±0.001 


16.012± 0.001 


18.345 ±0.002 


5546.026 


17.412± 0.001 


18.298 ±0.002 


16.088 ±0.001 


16.010± 0.001 


18.267 ±0.002 


5654.585 


17.242 ±0.001 


18.023 ±0.002 


16.107 ±0.001 


16.026 ±0.001 


18.397 ±0.002 


5729.072 


17.207 ±0.001 


18.258 ±0.002 


16.113±0.001 


16.048 ±0.001 


18.645 ±0.003 


5877.216 


17. 147 ±0.001 


18.047 ±0.002 


16.090 ±0.001 


16.013 ±0.001 


18.219 ±0.002 



NOTE. — Light curves are in ST magnitudes. In this filter, the offset from ST to AB magnitudes is m^B ■ 
msT = 1.532. 



at a time rather than pairwise, and second, we add in quadra- 
ture an extra uncertainty of 0.03 mags to each image. This is 
about the level of the high-frequency variability seen in image 
A in these data, and interpreted variously as microlensing of 
relativistic knots (Schechter et al. 2003), echoes of the intrin- 
sic quasar variability due to luminosity-dependent accretion 
disk area (Blackburne & Kochanek 2010 1, or microlensing o f 
stochastic hot spots on the disk (e.g., Dexter & Agol 2011). 
We do not currently have the capability to test these hypothe- 
ses (at least, not while running a multiwavelength simulation), 
so we add the extra uncertainty to avoid skewing our results 
by trying to fit this variability with stellar microlensing. 

Our B-band data includes flux from the redshifted Lya 
emission line, and R, H, and K are also somewhat affected 
by the Cm], H/3, and Ha lines, respectively. But in these 
broad filters, the line flux only amounts to a few percent of 
the total flux, and we do not expect it to ha ve a strong effect 
on the microlensing results (D ai et al.|2010) >. 

The light curve of the magnitude difference of the two 
quasar images, with error bars as described in this section, 
is shown for each band in Figure [T] together with an example 
of one of the many microlensing light curves that we fit to the 
data. Since this is the flux ratio of the two quasar images, and 
is corrected for the time delay, the variability is due to mi- 
crolensing. The trend of increasing microlensing variability 
with decreasing wavelength, which is qualitatively expected 
for an accretion disk, can be easily seen. 

2.2. UVData 

We have observed HE 1 104 in the F275W filter (rest-frame 
0.083 /im) using the UVIS channel of the Wide Field Camera 
3 aboard HST at 10 epochs roughly evenly spaced between 
2009 December 4 and 2011 November 11. At each epoch, 
we combine four images totalling 2524 seconds of exposure 
time using the Multidrizzle task within the PyRAF software 
packagq^j We take care to prevent the automatic cosmic ray 
rejection from masking out the PSF cores of the quasar com- 
ponents or bright comparison stars, manually removing some 
pixels from the bad pixel mask before repeating the final driz- 
zle. The resulting images are clean and contain little besides 
the two quasar images and a handful of other pointlike ob- 
jects. 

We use aperture photometry to determine the fluxes of the 
quasar images and three other comparison objects. This is 
a very simple process because of the paucity of sources and 

6 PyRAF and Multidrizzle are products of the Space Telescope Science 
Institute, which is operated by AURA for NASA. 



relatively large separation between the quasar images. We 
use a square l."6 aperture, and measure the sky background 
in a square annulus just outside the aperture, with an outer 
diameter twice that of the aperture. We then use the header 
keywords PHOTFLAM and PHOTZPT to convert the instru- 
mental magnitudes to the ST system (in this filter, the offset 
from ST to AB magnitudes is mAB-^ST = 1 -532). Tablepllists 
the ST magnitudes and formal uncertainties of quasar images 
A and B, as well as three comparison objects that we label 
a, b, and c. The first two of these, a and b, are the bright 
stars northeast of the lens labeled 4 and 3 (respectively) by 
Wisotzki et al. ( 1995 i, and the last is located at a position 
(Aq cos d , Ad) = (-55.7, -38 ."2) relative to quasar image A. 
Objects a and b show little variability, while object c varies by 
~0.3 mags and may well be a quasar. 

Since HE 1 104 has a relatively high redshift, these UV ob- 
servations probe the far-UV region of the quasar spectrum be- 
tween the Lyman limit and the Ly/3 emission line, from A ~ 76 
to 91 A in the rest frame. Though there may be some contami- 
nation from the remainder of the Lyman series, this is unlikely 
to be strong, and we consider the majority of the UV flux to 
be coming from the quasar accretion disk. 

We shift the light curve of image B and use linear interpo- 
lation to estimate its delay-corrected light curve, resulting in 
a 7-epoch light curve. We add 0.05 mags in quadrature to the 
uncertainties of each image to account for the systematics in- 
troduced by this process. Fortunately, the mismatch between 
each measurement of image A and the nearest shifted mea- 
surement of B is usually fairly small, on the order of 10-30 
days. The UV flux ratios are shown in the top panel of Fig- 
ured] 

2.3. X-Ray data 

For our investigation of the X-ray properties of the quasar, 
we use the X-ray light curves of HE 1104 from Chen et al. 
( |2012b| l and |Chartas et al.| ( |2009l >. These data consist of soft- 
band (0.4-1. 2 keV) and hard-band (1.2-8 keV) absorption- 
corrected count rates at nine epochs. We convert the count 
rates to instrumental magnitudes and symmetrize the result- 
ing error bars by taking the geometric mean of the upper and 
lower errors. 

The sampling of the X-ray light curves is generally not 
well-matched to the time delay between the two quasar im- 
ages. In particular, the first epoch is separated from the rest 
by a period much longer than the delay. For this epoch we do 
not apply any delay correction, but simply augment the error 
bars in quadrature by 0.078 mags as previously described. For 
the remaining epochs we linearly interpolate the light curve 
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FIG. 1 . — Time delay-corrected difference in magnitudes between the A and B images of HE 1 104 in each filter as a function of time (points with error bars). 
Some systematic errors have been added as described in Section[2] Curves show one of the best-fitting out of millions of microlensing light curves found by our 
simulation. The soft X-ray light curve is in its own panel because of its much larger range. The hard X-ray data are nearly identical, and we do not show them. 



of image B, then add the same extra uncertainty, except for 
the last two epochs. For these epochs, the observational ca- 
dence is within 10 days of the time delay, so that we can easily 
match one observation of B with the next observation of A. 
We add only half the extra uncertainty (0.039 mags) in these 
cases. The resulting soft-band X-ray flux ratios are shown in 
the bottom panel of Figure [1] 

Unlike some o ther lensed quasars s uch as Q 2237+0305 
HChen et al.|20Tfl > and RXJ1 131-1231 ( [Chartas et allpOH) , 
HE 11 04 does not show significant differences between its 
hard and soft X-ray light curves. This immediately indicates 
that our microlensing simulation analysis will not be able to 
give definitive evidence of a difference in source size for these 
two bands. However, we ought to be able to place upper limits 
on the magnitude of the logarithm of their size ratio, which is 
an interesting result in its own right. 

3. MICROLENSING SIMULATIONS 

We model the gravitational microlensing of the two quasar 
images using the Bayesian Monte Carlo technique described 
by Kochanek (2004) and updated by Poindexter & Kochanek 
( 2010a|bp . This method uses ray-tracing to generate magni- 
fication patterns that encode the microlensing magnification 
experienced by each quasar image on top of its magnifica- 
tion from the "macro" -lensing by the lens galaxy as a whole. 
We generate new patterns for each epoch in the quasar's light 



curve, allowing the stars to move between epochs. This 
causes the pattern of high-magnification caustics and low- 
magnification troughs to evolve. Together with the bulk rel- 
ative motion of the quasar and the lens galaxy, this leads to 
variability in the microlensing magnifications. This work is 
the first to simultaneously fit observations at a variety of wave- 
lengths and include the effects of stellar motions. 

To generate the magnification patterns, we scatter stars ran- 
domly across a portion of the image plane. The mass func- 
tion that we use for the microlens stars is a power law with a 
slope dlnN/dlnM = —1.3 and a mass range M max /M m j n = 50. 
We vary the mean mass (M) between 0.03M Q and 10M Q 
in six logarithmically spaced steps. At the location of each 
quasar image, the fraction of the surface mass density made 
up of stars (as opposed to smoothly distributed dark mat- 
ter) is determined from a global parameter, fujL- This pa- 
rameter varies between (no stars) and 1 (all stars), and 
specifies the amplitude of the stellar component (in a lens- 
ing model consisting of a de Vaucouleurs stellar component 
and a Navarro-Frenk- White dark matter component) relative 
to the best-fitting mo del with only stars. These models come 
from the work of Poindex ter et al.| ( |2007| >, and are constrained 
not only by the relative positions of the quasar images and 
lens galaxy, but by the images' near-IR fluxes, which are not 
strongly affected by microlensing, and by their estimate of 
the lensing time delay. They find that a value of f M n =0.3 
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is most likely; this estimate would probably not change much 
if more recent time delay estimates were used instead. In the 
interest of fair sampling, we allow f M / L to take the values 0.1, 
0.3, and 1.0. From these same lens models we obtain the total 
lensing convergence n tol and shear 7 at the positions of the 
quasar images, which together with (M) and f M / L determine 
the overall character of the magnification patterns. We set the 
rest-frame one-dimensional velocity dispersion of the stars to 
301 kms" 1 ; this is calculated from the monopole strength of 
a singular isothermal sphere plus external shear (SIS7) model 
for the lens potential. This has been shown to be a good esti- 
mator for the stellar velocity dispersion (Treu et al. 2006). 

The effect of the quasar's finite size is simulated by 
convolving the magnification patterns with a source light 
profile. We use a profile with surface brightness I(r) oc 
[exp[(r/rA) 3//4 ] - 1] _1 , where r\ is defined by the relation 
kT e ff{r\) = hpc/X. In standard thin disk theory, 



r\ = 



_ /45GA 4 M BH M 



V I6ir 6 h p c 2 



1/3 



4/3 



M, 



BH 



(9.7xl0 15 cm)f— ^1 ( „ 



2/3 



r]L E 



1/3 



(1) 



where A is the rest wavelength, L E( jd is the Eddington lumi- 
nosity, and r\ is the radiative efficiency. In this model, the 
disk is a multicolor blackbody with temperature T e ff oc r - ' 3 
and /3 = 3/4. We neglect the effects of the inner edge of the 
disk. We vary the projected area of the disk at a fixed wave- 
length corresponding to the observed-frame R band, allowing 
it to take values separated by 0.2 dex. We present most of 
our results, however, in terms of the half-light radius of the 
disk r\ji, which is proportional to the square root of the area 
divided by the cosine of the inclination angle, and is equal to 
2.44r\ when j3 = 3/4. The half-light radius is assumed to vary 
as a power law with wavelength, r\ji oc A^, and we allow its 
power-law slope £ to vary as well. Generically, £ is the re- 
ciprocal of the temperature slope, so for a thin disk its value 
is 4/3. We allow it to take values ranging from -1 to +2 in 
7 steps. Allowing the slope to vary is inconsistent with our 
choice of j3 = 3/4 for the light profile of the source, but fortu- 
nately the details of the light profile, apart from the half-light 
radius, do n ot have a strong effect on the microlensing (Mor- 
|tonson et al.|2005) . Since the X-ray flux does not arise from 
the same blackbody source as the UV/optical flux, we do not 
assume that the X-ray size lies on the same power-law curve, 
instead allowing the area of the source in soft and hard X-rays 
to vary independently. We do this by pursuing an iterative 
strategy, first performing an initial simulation that is ignorant 
of the X-ray data, and then following up with secondary sim- 
ulations for soft and hard X-ray data, each of which makes 
use of the output of the preceding simulation. The secondary 
simulations do not generate random trajectories for the quasar 
across the magnification patterns; instead they use the trajec- 
tories that successfully fit the data in the previous simulation, 
and they only vary a single new parameter: the area of the 
X-ray source. We do not take into account the lensing of the 
X-ray flux by the black hole, though given the compactness 
of the X-ray source it may be advisable to include this effect 
in future work ( |Chen et al.|2012a| ). 

We vary the inclination of the accretion disk between cos i = 
1 (face-on) and cos/ = 0.2 (nearly edge-on) in 5 steps evenly 



distributed in cos;; this is equivalent to a uniform prior on 
the orientation of the disk in three dimensions. Since an in- 
clined disk looks like an ellipse to the observer, we also let 
the position angle of the projected disk take 9 values evenly 
distributed between <\> a = 0° (major axis aligned with north) 
and cf) a = 160° (rotated east of north). These parameters have 
subtle effects on the smoothed magnification patterns, and 



Poindexter & Kochanek (2010a) show that with moving pat- 
terns it is possible to constrain their values. But we note 
that they are helped by the relative importance of the ran- 
dom stellar motions in the lens system they consider (Q 2237), 
and the only other attempt to use this method so far does 
not obtain stron g constraints on these parameters for the lens 
HE 0435-1223 ( |Blackburne et al.|2011a| > 



As described by Kochanek (2004) and Poindexter & 
Kochanek (2010a b), this analysis method generates billions 
of trial paths that the quasar may take across the magnifica- 
tion patterns, and evaluates the likelihood of each trial based 
on the goodness of fit of the resulting light curve. The calcu- 



lations are d escribed in Equations (4) and (5) of Blackburne 
(201 la I; we set the "rescale" parameter f£ to 2.0. 
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1.0, 

and 1 .0 tor the UV/optical, soft X-ray, and hard X-ray bands, 
respectively. Since each trial is associated with some loca- 
tion in parameter space, this Monte Carlo sampling gives us a 
reasonable idea of the joint likelihood function of our param- 
eters. We then multiply the likelihood by our priors to obtain 
a Bayesian posterior probability density functions (PDFs) for 
our parameters. Our priors are mostly uniform, either in lin- 
ear or logarithmic space. The mean microlens mass (M) and 
the variable parameterizing the fraction of the surface den- 
sity in stars f M i L have logarithmic priors, and the inclination 
cos/, the disk position angle <p a , and the wavelength slope 
£ = d\nri/2/dh\ A have linear priors. For the half-light radius 
r\j2 itself we use both linear and logarithmic priors, because 
they can both be s ensibly applied (see discussion by Black- 
|burne et al.|201 la) , and in order to evaluate the robustness of 
our results to changes in this prior. For the starting positions 
of the quasar on the magnification patterns we choose a simple 
uniform prior, and for its velocity we choose a circular Gaus- 
sian. The velocity prior is calculated in the manner described 
by Blackburne et al.| l |2011a) , centering the Gaussian on the 
projection onto the lens plane of the velocity of the solar sys- 
tem relative to the cosmic microwave background (CMB), and 
setting its width to the quadrature sum of estimates of the pe- 
culiar velocities of the lens galaxy and the source galaxy. All 
velocities are corrected for cosmological time-dilation, and 
are converted to angular velocities using the angular diame- 
ter distances Dol, ^os, and Z\s (observer to lens, observer 
to source, and lens to source, respectively). Projected into the 
lens plane, the effective one-dimensional width of the velocity 
prior is 

l+z L \ 2 (Do h \ 2 

l+z s ) \D 0S ) 



J eff 



Or +0\ 



(2) 



We us e the prescription cited by Mosquera & Kochanek 
( |2011[ ) for the estimates of the peculiar velocity dispersion, 
and find ctl = 275 km s" 1 and erj = 201 km s" 1 , yielding a width 
for the velocity prior of a e s = 290 km s" 1 . This value is also 
approximately the typical speed of the quasar across the pat- 
tern. 

4. RESULTS 

The product of the likelihood distribution with our priors 
gives a joint posterior probability distribution for our parame- 
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FIG. 2. — Filled contours indicate the 68%, 95%, and 99% confidence levels 
of the velocity prior described in SectionR] Open contours show the poste- 
rior probability distribution of the transverse lens velocity (relative to that of 
the source and the observer). The direction of motion of the lens galaxy is 
120 ± 40 degrees, with a 180-degree degeneracy due to the symmetry of the 
elliptical disk model. 

ters. In this section, we examine the projections of this distri- 
bution, specifically the posterior PDFs for the velocity of the 
quasar relative to the lens galaxy, the inclination of the accre- 
tion disk and the position angle of its major axis, the mean 
mass of the stars causing the microlensing, the half-light ra- 
dius of the quasar in the optical, soft X-ray, and hard X-ray 
bands, and the power-law slope of the half-light radius of the 
accretion disk with wavelength. 

4.1. Velocity 

The moving magnification patterns break the strict degener- 
acy between the velocity of the quasar relative to the lens and 
the mean mass of the stars, and allows us to put some actual 
constraints on the velocity. The posterior PDF for the trans- 
verse velocity of the lens (relative to the observer and source) 
is shown, together with our prior on the same quantity, in Fig- 
ure [2] More precisely, this quantity is the angular velocity of 
the source across the magnification pattern, multiplied by the 
lens distance Dol and time-dilated to the lens redshift, with a 
change of sign to make it a lens velocity rather than a source 
velocity. This differs somewhat in magnitude from the trans- 
verse peculiar velocity of the lens galaxy because it includes 
terms from the source and observer motions, but the lens term 
dominates and the difference is therefore small. Though the 
radial profile of the probability distribution is mostly deter- 
mined by the prior, we do see a preference for a direction of 
motion of roughly 120 ±40 degrees (68% confidence), with 
a 180-degree degeneracy related to the bilateral symmetry of 
the (generally inclined) accretion disk. 

4.2. Disk Inclination 

The posterior PDFs for the accretion disk inclination cos i 
and the position angle of its major axis <\> a are shown in Fig- 
ures[3]and|4] A face-on disk is favored, with cos / = 1 .0 about 
four times as likely as cos/ = 0.5. If we assume that the disk 




0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 
cos; 

FIG. 3. — Posterior distribution for the accretion disk inclination cosi. A 
face-on disk has cos i = 1 . 
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FIG. 4. — Posterior distribution for the accretion disk major axis position 
angle </>„, measured in degrees East of North. Face-on (cos; =1) solutions 
are excluded in calculating this result. 

is in fact inclined and discard trials with cos i = 1 .0, then its 
major-axis position angle is fairly well-determined, with a 
value of 150 ±33 degrees East of North (68% confidence). 

4.3. Mean Microlens Mass 

FigureBJshows the posterior PDF for the mean mass (M) of 
the stellar microlenses. This quantity determines the scaling 
between the natural angular units of the magnification pat- 
terns (microlens Einstein radii) and units of physical length 
(e.g., cm) in the source plane. Even with our moving magni- 
fication patterns, it is still somewhat dependent on the veloc- 
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ity prior, in the sense that the very large transverse velocities 
suppressed by the Gaussian wings of our velocity prior cor- 
respond to very large mean masses. But we are confident in 
the appropriateness of this prior; in any case it seems quite 
unlikely that the mass function of stars in the lensing galaxy 
has a mean higher than 1 M©. The distribution is fairly broad, 
and favors mean masses of 0. 1 to 0.3 M©. 

4.4. Accretion Disk Size 

The posterior PDFs for the deprojected half-light radius 
of the accretion disk in the R band, and of the soft (0.4 to 
1.2 keV) and hard (1.2 to 8keV) X-ray emission regions are 
shown in Figure [6] The logarithm of the /?-band radius (mea- 
sured in cm) is lo.O^gj with the logarithmic prior. With the 
linear prior, this value rises to 16.2 ±0.3. The X-ray PDFs do 
not converge at small sizes, as the finite resolution of the mag- 
nification patterns prevents us from going to smaller sizes, but 
the cutoff scale is smaller than the gravitational radius of the 
black hole (regardless of which black hole mass estimate is 
chosen). The three vertical lines in Figure [6] show the gravi- 
tational radius, the /?-band radius estimated from the flux of 
the quasar, and the /?-band radius predicted by the thin disk 
model, assuming an accretion efficiency r\ = 0. 1 and an Ed- 
di ngton fraction L/L sdd = 1 ■ All assume the H/3 mass estimate 
of Assef et al. ( 2011| l. With the logarithmic prior the source is 
smaller than log(ri/2/cm) = 15.33 at 95% confidence in both 
soft and hard X-rays. With the linear pri or, the upper limits 
are 15.71 (soft) and 15.59 (hard). If the |Peng et"aL| (2006 ) 
black hole mass estimate is correct, then half the X-ray light 
is originating within a radius less than about 6 (logarithmic 
prior) or 12-14 (linear prior) gravitational radii. 

The deprojected half-light radii are calculated from the pro- 
jected area of the source A = 7r(r 1 /2/2.44) 2 cos/. With the 
logarithmic (linear) prior, we find that log(A/cm 2 ) = 31.54!]j" 
(32.26^5o). Likewise, the 95% confidence upper limit on the 
logarithm of the source area is 30.23 (31.34) for soft X-rays, 
and 30.20 (31.11) for hard X-rays. These projected areas 
are more appropriate for comparison with previous work that 
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FIG. 6. — Posterior distribution for the half-light radius of the quasar accre- 
tion disk in the R band and the two X-ray bands. In the upper panel we use 
a logarithmic prior, and in the lower panel we use a linear prior. The stellar 
mean mass (M/Mq) is fixed at 0.3. The gravitational radius of the black 
hole and the R-band size estimates based on the quasar flux and the thin disk 
model are shown as a vertical lines, marked "Grav," "Flux," and "Disk" re- 
spectively. The "Grav" and "D isk" values assume the H/3 black hole mass 
estimate of Assef et al.|J2011) , and can shift up to ~0.5 dex based on the 
uncertainties or by using the Ha or Civ estimates. 

does not vary the inclination of the source, but since HE 1 104 
seems to have a relatively face-on disk, this correction is not 
very important. 

The soft and hard X-ray half-light radii PDFs are nearly 
identical, which is consistent with our expectations given the 
similarity of their light curves. To investigate the question of 
which is larger, we show in FigurelTlthe posterior PDFs for the 
logarithm of the ratios between the sizes in the various bands. 
The solid curve indicates the ratio of the hard band size to the 
soft band size. It peaks fairly sharply near zero, indicating a 
ratio of unity. The symmetric tails on either side mean that we 
cannot with our current data distinguish which is the larger, 
but we can say that their sizes do not differ too much: at 68% 
confidence the logarithm of the ratio falls between -0.45 and 
+0.46. Figure I7J also indicates that the X-ray sizes are much 
smaller than the 7?-band sizes, with Alog(ri/2) = -1 .40!J;« for 
the soft band and— 1.35!SJ in the hard band (68% confidence). 
These constraints are tighter than the previous figure would 
imply, due to covariances between the optical and X-ray sizes 
that are not apparent in Figure [6] 

Figure [8] shows the posteriorPDF for the power-law slope 
of the wavelength dependence of the half-light radius, a 
quantity we call £. For an accretion disk that radiates as 
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FIG. 7. — Posterior distributions for the ratios of the quasar's half-light ra- 
dius. The black curves indicate ratios of the X-ray sizes; for these curves the 
mean mass (M) is set to O.3M0. The X-ray emission is much more com- 
pact than the observed-frame ft-band emission (labeled "Optical"). Also, our 
X-ray data rule out a very large difference in the sizes of the hard and soft 
X-ray sources. The gray curves indicate the size of the accretion disk at the 
(observed-frame) UV and /-band wavelengths, relative to the S-band size. 
Since we parameterize the disk using the S-band size and the wavelength 
slope f, the gray curves are simply scaled versions of the posterior distribu- 
tion for £ (see Figurepj). 
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FIG. 8. — Posterior probability distribution for the power-law slope of the 
wavelength dependence of the half-light radius. The standard thin disk model 
predicts a value of 4/3 since T e g <x r~ 3 ' 4 . 

a multi-temperature blackbody, this exponent is the recip- 
rocal of the power-law slope (3 of the temperature profile 
r e ff(r) oc r"* 3 . So for the standard thin disk model, we ex- 
pect £ = of log r^/^ log A = /3" 1 = 4/3 (neglecting the effect of 
the inner disk edge, as is reasonable to do at radii of many 
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FIG. 9. — Joint posterior probability distribution for the mean microlens 
mass (M) and the projected source area A in the R band (top panel), soft X- 
rays (middle panel), and hard X-rays (bottom panel). The solid black (dashed 
gray) curves indicate a logarithmic (linear) prior on the area. The vertical line 
indicates the square of the black hole's grav itational radiu s, as calculated us- 
ing the black hole mass estimate of Assef et al. 12011). The gray region 
cannot be sampled because of the finite resolution ot the magnification pat- 
terns. 

r g ). Our probability distribution peaks at a value of 1.0, with 
a median value of 0.84, and with 68% of the probability lying 
between values of 0.44 and 1 .30. This is consistent with the 
thin disk value, though it is interesting that the majority of the 
probability lies toward smaller values of £, or (equivalently) 
steeper temperature profiles. Using this posterior distribution, 
we show in Figure [7] the logarithm of the size ratio between 
the observed-frame UV and R-band disks, and between the J- 
band and /?-band disks. These are just scaled versions of the 
distribution in Figure [8] 
Figures|^and|7]andtne source size estimates quoted in this 



THE STRUCTURE OF HE 1 104-1805 FROM INFRARED TO X-RAY 






2.0 



1.5 



S- i.o 






a 
o 
a. 
x 
W 



0.5 



0.0 



-0.5 



1 


1 1 ^ 




i 










/ 








/ 




A s 










r 1 \ 


\ \ 




/ 








\ \ ' 


- 




/ (' 




/ 


V \ ' 






\ I ' 


n 


"~~\ 


\\ V i 






\ ' / 
/ / ' / 
/ / ' / 
/ / ' / 


i 
i 
i 


/ 
\ 
\ / 


\ \ \ y i 

/ / / // / 

i in 

Jr ' 






( \ v v. 




\ J 


f 1 i 
I i 






\ V^ x 




\*^s 


-v ' 




ll 




— ^ 








i i 


i 


1 



28 



29 30 31 32 33 
Projected Area log(A/cm 2 ) 



34 



FIG. 10. — Joint posterior probability distribution for the exponent £ and 
the projected R-band source area A. The solid black (dashed gray) curves 
indicate a logarithmic (linear) prior on the area. The gray area is excluded 
because the observed-frame UV disk would be smaller than the resolution of 
the magnification patterns. 

section are produced with the mean microlens mass (M) fixed 
at its most likely value of 0.3M Q . We do this rather than 
marginalizing over (M) in order to avoid an artificial decrease 
in the probability of small source sizes due to the finite res- 
olution of the magnification patterns. This artificial decrease 
happens because there is no contribution to the probability of 
small sizes from trials with large (M) (and thus with a pro- 
jected pixel size larger than the area of the source). To ex- 
plore the effect of this choice, we show in Figure [9] the joint 
posterior PDF for the mean mass (M) and the projected source 
area A = ir(r l / 2 /2A4) 2 cosi. Some covariance between these 
parameters can be seen; this is expected because the area is 
measured in units of the square of the microlens Einstien ra- 
dius. The figure also shows the resolution limit of the magni- 
fication patterns, and it is clear that the /?-band source is well- 
resolved in all cases, but that the X-ray sources are not much 
larger than the pixels in the magnification patterns (despite the 
higher-resolution patterns we use in these simulations). 

Another concern for the accretion disk simulations is that 
at large values of £ the HST UV observations (i.e., the bluest 
apart from X-rays) will probe source sizes small enough to be 
unresolved, even though the source is resolved in the R band. 
Figure 10 shows the joint posterior PDF of the exponent £ and 
the /?-band projected source area A. We also show the region 
where the UV source is unresolved, with (M) = 0.3Mq. The 
probability distribution clearly converges at small areas, so 
this concern is unfounded. 

5. DISCUSSION AND CONCLUSIONS 

HE 1104 is the third lensed quasar to be analyzed using 
dynamic microlensing magnification p atterns, after Q2 237 
(Poindexter & Kochanek 2010a b; Mosquer a et al.|2013) > and 
HE0435-1223 (Blackbur ne et al.||2011a| ). The inclusion of 
the random motions of the stars allows us to constrain the 
inclination of the accretion disk, the position angle of its pro- 
jected major axis, and the direction of the lens galaxy's mo- 



tion relative to the quasar. In this case, the data favor a low 
disk inclination, with cos i = 1.0 about four times as likely 
as cos; = 0.5 (se e Figure [3]). This is similar to the result that 
Poindexter & Kochanek ( 2010a) find for Q 2237, and supports 
the "unification" model for AGNs, which predicts that bright 
Type 1 quasars such as these have a low inclination. For mod- 
els with a nonzero inclination, we find that the major axis of 
HE1104's projected disk is about 150 ±33 degrees east of 
north (see Figure |4] We also find that the lens galaxy's mo- 
tion is likely toward the northwest or southeast, with a most 
likely angle of 120 degrees east of north (or 30 degrees west 
of north), and with 68% confidence error bars of 40 degrees. 
The magnitude of the velocity is unfortunately not robustly 
constrained, and depends on our velocity prior (see Sectionp] 
and Figure [2]). 

We also produce a posterior PDF for the mean mass of the 
stars in the lensing galaxy (see Figure K). The distribution 
peaks between 0.1 and 0.3M Q , which is consistent with what 



is observed in other lens gala xies (Poindex ter & Kochanek 
2010bHBlackburne et al.|201 la) and in the Milky Way (|HoItz- 



man et al. 1998; Zoccali et al. 2000). Our distribution implies 
a slightly higher mean mass than that of Chartas et al. (2009), 
which is probably due to differences in the velocity priors. 
Their analysis uses static magnification patterns, and (prob- 
ably more importantly) their peculiar velocities are smaller 
than the ones that we have chosen, which would tend to make 
their mean mass estimates smaller. 

Our results for the half-light radius of the accretion disk 
in the observer-frame R band (0.2 //m in the rest frame) 
agree with t hose of p revious studi es of HE 1 104 
etaLl[2008l |Morgan et af. 2008 [ Chartas et al.l |2009 ) 
these studies, we find that HE 1 104, like several other lensed 
quasars, has an accretion disk much larger than would be ex- 
pected from either the thi n disk model or the q uasar flux (see 
Equations 2 and 3 of Poindexter et al.| ( |2008) ). These es- 
timates, marked "Flux and "Disk" respect ively, are plotted 
for the observed-frame R band in Figure [6] Poindexter et al. 



(2008) suggest a small value for the accretion disk tempera- 



ture profile slope /3 as a possible solution to the discrepancy 
between the microlensing size and the flux size, supported by 
their estimated value of /3 = 0.61!Jn, smaller than the canon- 
ical 0.75. But our result for the slope of the size-wavelength 
relation, £ = /T 1 = LOS, implies a rang e of 0.77 to 1.85 
for j3. This is roughly consistent with the Poindexter et al. 
(2008) result, but some tension remains. It seems that our 
new data indicate that the area of the disk changes less with 
wavelength than was previously thought. Indeed, our value fa- 
vors a temperature profile steeper than the standard thin disk 
model. This would only serve to ex acerbate the flux size dis- 
crepancy (see Morgan et al.||2010 l. Multiwavelength obser- 
vations of more lensed quasars will shed more light on this 
ongoing mystery. 

Our upper limits on the size of quasar at X-ray wavelengths 
are quite strong. It is clear that the majority of the hard and 
soft X-ray flux is coming from the innermost ^ 30 r s , assum- 
ing the H/3-based black hole mass estimate of |Assef et al. 



(201 1 ). This result is comparable to the upper limit that|Char- 



tas et al. (2009 ) find using a single X-ray data point. This lack 
of improvement is partly due to our practice of dividing the 
X-ray flux into two bands and of marginalizing over a variety 
of source inclinations (both of which broaden the posterior 
PDF), and partly due to the fact that we have not yet been 
able to sample the X-ray light curve on time scales smaller 
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than the source crossing time (see Mosquera & Kochanek 
201 1] >. For HE 1104 the soft and hard X-ray data are very 
similar, and so we are unable to determine which band has 
a larger size, though we do conclude that their sizes do not 
differ by more than 0.46 dex. This result is similar to those 



of Blackburne et all (12011 all for HE 0435 and Morgan et al. 

,012) 



( |2012h for Q JO 158-4325. The question of the two X-ray 
bands relative sizes is quite interesting, since it can address 
the corona/reflection model. If the direct component of the 
emission dominates, the hard band, coming from a hotter elec- 
tron population, ought to be smaller. But a prominent reflec- 
tion from the accretion disk could reverse this result, since the 
reflected spectrum is harder than the input spectrum, and may 
well be more extended. Recent papers presenting X-ray ob- 
servations have used simple microlensing models to address 



this question ( jChartas et al. 2012 Chen et al. 2012b I, but 
full microlensing simulations are needed to put rigorous con- 
straints on the sizes. Some of the X-ray data display fairly 
strong chromatic variation, so interesting results should be 
forthcoming. 
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